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Experimental Measurements of Atmospheric Neutrinos 

Edward Kearns a * 

a Department of Physics, Boston University 
590 Commonwealth Ave. 
Boston, MA 02215 

This talk reports the latest indications of an anomaly in the measurements of atmospheric neutrinos. New 
results from Soudan-2 and Super-Kamiokande provide evidence that the ratio of to z/ e interactions is not as 
expected. High energy Super-Kamiokande data indicates the cause is a deficit of upward-going and the zenith 
angle dependence of the effect is consistent with neutrino oscillations. Upward-going muon measurements by 
several detectors are discussed, but in total they provide inconclusive evidence for the anomaly. 



1. Introduction 

Large underground detectors originally built to 
search for proton decay are also exposed to a flux 
of neutrinos created by cosmic ray showers in the 
upper atmosphere. Neutrino interactions in the 
detector can mimic proton decay, and therefore a 
great effort has gone into predicting the rate and 
topology of the neutrino background. A byprod- 
uct of this effort has been the recognition of an 
anomaly: the relative rate of and v e interac- 
tions disagrees with expectation, and the baseline 
and energy dependence of the disagreement sug- 
gests that neutrino oscillation may be the culprit. 

The nature of these experiments is to measure 
the rate of neutrino interactions in the detector 
and to compare that rate with theoretical pre- 
diction. The theoretical task is to predict the 
neutrino flux as a function of energy, direction, 
and flavor, taking into account measurements of 
cosmic ray flux, geomagnetic cutoff, and produc- 
tion and decay of secondary mesons. Detailed 
Monte Carlo programs are used to estimate the 
cross section for neutrino interaction and simulate 
the response of the detector. The experimental 
task is to measure as much as possible about the 
neutrino interaction, in particular, the energy, di- 
rection, and flavor of the final state lepton, from 
which the neutrino properties are inferred. 



2. Flux prediction 

One of the tenets of the atmospheric neutrino 
anomaly is that the flux raticQ (v^/ve) is more 
accurately predicted that the or v e flux alone. 
The principal effect is that cosmic ray showers 
consist mostly of pions, which decay to -t- z^, 
and the \x decays to to e + v e + i>^, resulting in 
a flux ratio (v^/v e ) ~ 2. The authors of detailed 
calculations of the flux models have collaborated 
to compare results [Q and reached an understand- 
ing of many of the differences between their ear- 
lier publications. Two updated calculations ||[| 
cover a wide energy range (10 MeV-10 TeV) and 
are used by current experiments. There are now 
also an assortment of new high altitude /i mea- 
surements available for comparison [|]-[| . Further 
details of flux calculations were presented at this 
meeting by T. Stanevj7|. However, there is cur- 
rently no indication that poor knowledge of the 
predicted flux could be responsible for the exper- 
imental anomalies described below. 



3. Summary of results 

To study the {v^/v e ) flux ratio, most ex- 
periments calculate the double-ratio or ratio-of- 
ratios: 
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♦Presented at TAUP97, the 5th Int'l Workshop on Top- 
ics in Astroparticle and Underground Physics, Sept. 7-11, 
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{NJN e ) MC 
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technically [y^ + ^ M )/(i/ e + v e ) 
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where N refers to the number of events where the 
final state lepton is classified as /u-like or e-like 
by some identification algorithm. As mentioned 
above, the ratio of f M to v e flux is accurately pre- 
dicted; in addition, other theoretical and experi- 
mental uncertainties largely cancel. Table |l| lists 
previous measurements^ of R for E v ~ 1 GeV, 
including the new results from Soudan and Super- 
Kamiokande discussed in this paper. The kine- 
matic limits differ somewhat from experiment to 
experiment, with minimum lepton momenta re- 
quirements from 100 to 200 MeV/c. Kamiokande 
restricted their sample to^| E V i S < 1.33 GeV, 1MB 
to p < 1.5 GeV; the other experiments did not 
specify an upper limit, but all results are domi- 
nated by ~ 1 GeV neutrinos. 



Table 1 

Summary of R measurements, E v ~ 1 GeV. 

Experiment kt-yr events R (data/MC) 
Super-K|j] 
IMBgl 
Kamr|§ 
Soudan-2 jl| 
Frejuspi 
NUSEXp[ 



Kamiokande also studied events with E V i s > 
1.33 GeV (multi-GeV) and included partially con- 
tained (PC) events where a track was detected 
exiting the inner detector [^5). They measured a 
low value for R of 0.57±^;^f ±0.07, but more inter- 
esting was the dependence of R on zenith angle. 
Neutrinos that travelled ~ 10 4 km from below 
showed a small value of R, but those that trav- 
elled ~ 10 km from above agreed with expecta- 
tion, suggesting an oscillation length somewhere 
in between. 
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3 The first uncertainty quoted is statistical, the second is 
systematic; this convention will be used throughout the 
paper. 

4 Visible energy (E v i s ) is defined as the energy of an 
electromagnetic shower that produces a given amount of 
Cherenkov light. 



Most of the 1MB multi-GeV exposure had a re- 
striction on the maximum number of PMT hits 
(to concentrate on proton decay); the restriction 
was eventually removed, so they made a separate 
analysis |l6fl of their last 2.1 kt-years for E V i s > 
.95 GeV and measured R = 1.40±g;^j ± 0.14 with 
no zenith dependence. The caveats are: lim- 
ited statistics of 72 events (some overlapping with 
Ref.||); no outer detector to help identify PC 
muons; coarser sampling (4% photon coverage), 
resulting in only 90% correct e//x identification. 

4. Soudan-2 



Until recently, atmospheric neutrino results 
seemed to be divided between water Cherenkov 
detector s |9| ,|l0| (anomalous) and iron calorime- 
ters |l3| , |l4| (as expected). The Soudan-2 collab- 
oration, which operates a fine-grained iron track- 
ing calorimeter in Minnesota, U.S.A., has recently 
published results |pd| which support the anomaly 
seen in the water Cherenkov experiments. 

At this meeting, T. Kafka has updated the re- 
sults from Soudan-2 to 3.2 kt-years |T^|. They 
measure 91 single-prong track events (mostly 
charged current (CC) v^) with p > 100 MeV/c, 
and 137 shower events (mostly CC-^ e ) with p > 
150 MeV/c. The interaction vertex is allowed as 
close as 20 cm from the edge of the detector; with 
only 32 gm/cm 2 of shielding they observe a sig- 
nificant (25-30%) background from gamma rays 
and neutrons associated with nearby cosmic rays. 
However, they use an active shield of proportional 
tube planes lining the detector hall to veto most 
nearby cosmic rays, as well as separately estimate 
the remaining background rate as a function of 
flavor, depth into the detector, and energy. 

After background subtraction, the number of 
shower events matches their Monte Carlo predic- 
tion; however they observe 37% fewer tracks than 
they predict. Since the total flux is uncertain, it 
is better to consider the double-ratio, which they 
measure to be R — 0.61 ± 0.15 ± 0.05. Regarding 
the atmospheric neutrino anomaly, this is a con- 
siderable new piece of information, as the system- 
atics are very different from the water Cherenkov 
detectors. Although there are no demonstrated 
nuclear effects that would change the ratio of 
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to v e cross sections |17||, it interesting that the 
anomaly has also been seen in Fe as well as in 
H2O. What the difference is between these re- 
sults and those of Frejus and NUSEX (beyond 
what may be encompassed by large uncertainties) 
remains to be explained. 

5. Super-Kamiokande 

Super-Kamiokande is the next generation wa- 
ter Cherenkov experiment after 1MB and Kamio- 
kande. The detector resides nearby the old Kam- 
iokande detector in a mine near Kamioka, Japan. 
However, it is much larger (22.5 kton fiducial 
mass, versus 1 kton for Kamiokande and 3.3 kton 
for 1MB). It is instrumented with 11,146 PMTs, 
each 50 cm across, such that 40% of the inner 
surface area is active photocathode. The PMTs 
and electronics are of advanced design, with 2.5 
ns RMS timing for single photoelectrons. The in- 
ner detector is surrounded by an outer volume of 
water ~2.7 m thick that shields against incom- 
ing radioactivy and is instrumented with PMTs 
to tag penetrating muons. Further description of 
the detector, as well as new measurements of the 
8 B solar neutrino flux were presented at this con- 
ference by K. Inoue[p^[. 

The measurement from Super-Kamiokande is 
the result of a 414.2 live-day exposure (25.5 
kton-yrs) during the period from May 1996 to 
October 1997. The data is reduced from approx- 
imately 800K events per day to about 30 events 
per day by a series of software cuts. The most 
powerful requirement is the absence of hits in the 
outer detector, which indicates a fully contained 
(FC) interaction. The remaining events are then 
filtered by a visual scan, where the principal back- 
grounds are: (1) cosmic ray muons that evade the 
outer detector veto, typically by entering along 
cable bundles and then stopping in the detector, 
and (2) "flashing" PMTs that emit light due to in- 
ternal corona discharge. The partially contained 
sample is formed by a different reduction program 
from the same exposure^], since outer detector 

5 The results quoted here are updated from previous con- 
ference presentations, where PC livetimc was somewhat 
less than FC livetime, and PC data was scaled (by ~ 1.1) 
when FC+PC results were plotted. 



hits are now expected, and the background rejec- 
tion of entering cosmic rays is different. A 10.0 
live-year Monte Carlo sample of v interactions is 
passed through the same reduction chains, except 
for the visual scan. 

The remaining events, both data and Monte 
Carlo, are passed through the same reconstruc- 
tion code to: (a) fit the vertex of the interaction, 
by residual PMT timing, (b) count the number 
of Cherenkov rings, (c) estimate the direction of 
each ring, (d) estimate the energy of each ring, (e) 
determine the particle type (/i-like,e-like) for each 
ring, and (f) count the number of /i-decay elec- 
trons that follow each event. Most of the analysis 
is then done with the sample of events in which 
the number of rings found in (b) is exactly one. 
In most cases, this is the final state lepton from 
a charged current neutrino interaction; the prin- 
cipal contamination is single pion production as- 
sociated with neutral current (NC) interactions. 

The absolute energy scale was determined to 
~ 2.5% accuracy using several calibration signa- 
tures: LINAC electrons, radioactive sources, n°s 
and cosmic ray muons. About 9 photoelectrons 
are recorded for 1 MeV of visible energy. Conver- 
sion to lepton momentum takes into account the 
Cherenkov cutoff for muons. 

Data samples are defined using the same kine- 
matic criteria as in the Kamiokande experiment: 
Ve > 100 MeV/c, Pf, > 200 MeV/c, and E ms < 
1.33 GeV for sub-GeV; E ms > 1.33 GeV for 
multi-GeV FC. The partially contained sample 
is specified by a vertex in the inner detector and 
correlated hits in the outer detector; the mini- 
mum visible energy required is ~ 350 MeV. Be- 
cause (E v ) is >> 1 GeV for the PC sample, these 
data are added to to the FC multi-GeV sample; 
the CC lepton is assumed to be a muon, and no 
single-ring is required. The fiducial sample is re- 
stricted to events with vertex 2 m from the PMT 
wall (22.5 kton). 

The event totals are listed in Table || along with 
the totals for the Monte Carlo samples, scaled to 
25.5 kton-yrs. These yield the following values of 
R = {NjN e ) DATA /{NjN e ) M c-- 

m h rvv / °- 610 -o:o28 ± 0.049 (Honda flux) 
sub-UeV I 609 +o.O29 ± Q 049 (Bartol flux) 
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Table 2 

Event summary for 25.5 kt-year sample of fully-contained atmospheric neutrinos in Super-Kamiokande. 
Monte Carlo breakdown uses Honda flux. 



Monte Carlo prediction Monte Carlo breakdown 





Data 


Bartolg 


Honda@ 


CC-zv 


CC-i/ e 


NC 


sub-GeV e-like 


983 


788.9 


812.2 


2% 


88% 


10% 


sub-GeV /i-likc 


900 


1185.4 


1218.3 


96% 


0.5% 


4% 


sub-GeV multi-ring 


696 


753.7 


759.2 


43% 


24% 


33% 


multi-GeV e-like 


218 


190.9 


182.7 


7% 


84% 


9% 


multi-GeV F.C. //-like 


176 


229.7 


229.0 


99% 


0.5% 


0.4% 


multi-GeV P.C. (//-like) 


230 


305.0 


287.7 


98% 


1.5% 


0.6% 


multi-GeV multi-ring 


398 


450.1 


433.6 


55% 


30% 


15% 



multi-GeV 



0.659+' 
0.665 



-0.053 
+0.059 
-0.053 



± 0.081 (Honda flux) 
± 0.082 (Bartol flux). 



For both the high and low energy samples there 
is a significant deviation of the ///e ratio from the 
expected value of 1. The leading contributions to 
the systematic uncertainty in R are: {v^/v e ) flux 
(5%), neutrino cross section (4.6% for sub-GeV 
and 5.8% for multi-GeV), and single-ring selec- 
tion (3% for sub-GeV and 6% for multi-GeV). 

It is informative to check the relative rate of 
//-decay associated with the event sample. The 
decay electrons are detected as time separated 
hits from the neutrino interaction; most come ei- 



ther from associated tt + 



e + , or directly 



from /i^ 



in GG-v^ interactions. The per- 



formance of e/fi identification algorithms is no 
longer in question [|l9| , but the measured //-decay 
fractions check that the associated pion produc- 
tion is reasonably modeled in the Monte Carlo. 
Table || shows that the expected fraction of //- 
decay agrees well with the prediction; the fraction 
of //-decay found in stopping cosmic ray muons 
verifies the efficiency of the reconstruction. 

The Super-Kamiokande group had two inde- 
pendent analysis efforts that were used to check 
each other and minimize the possibility that some 
mistake would be made. The data were sep- 
arated after electronics calibration of the PMT 
data to photoelectrons and nanoseconds. Other- 
wise, everything was coded independently, includ- 
ing event reduction and reconstruction, Monte 
Carlo generation, and estimation of the energy 
scale. Beyond the independent code, the major 
differences in the second analysis were: (a) data 



Table 3 

//-decay fractions. 

Percentage of events with > 1 //-decay 

Data Monte Carlo 



stopping CR // 

//-like 

e-like 



74.0 ± 0.3% 
67.6 ± 1.6% 
9.3 ± 0.9% 



72.9 ± 0.4% 
68.1 ±0.1% 
8.7 ±0.3% 



Percentage of events with > 2 //-decay 
Data M.C. 



stopping CR // 0.0 ± 0.0% 
//-like 



e-like 



2.9 ± 0.6% 
0.2 ±0.1% 



0.0 ±0.0% 
4.1 ±0.1% 
0.1 ±0.0% 



reduction involved no scanning, (b) single-ring se- 
lection was based on an algorithm that classified 
events as single-ring or multiple-ring without at- 
tempting to count the number of rings, (c) e/// 
identification was somewhat simpler and less effi- 
cient (97% vs >99%), (d) some details of vertex 
and direction reconstruction were different. 

Upon comparison, the independent analyses 
were in exceptional agreement. Of the sub-GeV 
events found in the fiducial volume by the second 
analysis, 99.9% were found in the data sample of 
the first analysis, with 89% in the fiducial volume, 
consistent with the vertex fit resolution. Single- 
ring classification was in agreement 90% of the 
time. Comparing common events in the fiducial 
single-ring sub-GeV sample, vertices agreed to 84 
cm RMS, direction agreed to 2.5°, momentum 
agreed to 0.5%, 97% of the events agreed in parti- 
cle identification. The value of R for the sub-GeV 
sample of the second analysis is: 0.65±0.03±0.05; 
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the difference in value from the first analysis is 
understood to be due to differences in analysis 
methods and within their systematic uncertain- 
ties. In sum, the independent analysis provides 
reassurance that the deviation of R from unity is 
not due to experimental mistakes. 

To consider that the anomalous R is due to 
neutrino oscillation, one looks for a path length or 
energy dependence of the effect. The probability 
of two-flavor neutrino oscillation from v to v is 
given by: 

Am 2 (eV 2 /c 2 )I(km) 



Super-Kamiokande Preliminary 



P, 



sin 2 29 sin 2 1.27- 



(2) 



E(GeV) 

where 9 and Am 2 = |m 2 — m 2 , | are fundamental 
parameters that govern the neutrino mixing, and 
L and E are the path length and energy of the 
neutrino. The final state lepton direction and en- 
ergy are correlated with the incoming neutrino; 
for the sub-GeV sample, the mean opening angle 
for v^-ix is 54°, for v e -e it is 62°; for the multi- 
GeV sample it improves to < 15°. 

The samples are divided into 5 cos bins where 
is the angle between the outgoing lepton direc- 
tion and the nadirQ so down-going neutrinos that 
are produced directly overhead, with short travel 
distance, populate the bin near cos 6 = 1. Cal- 
culating R for each zenith bin results in Fig. |l|. A 
slight asymmetry is evident in the sub-GeV sam- 
ple, and a strong asymmetry is evident in the 
multi-GeV sample. If there were no anomaly, 
the R values would be around 1; for hypothet- 
ical oscillation parameters of sin 2 29 — 1 and 
Am 2 = .005 eV 2 /c 2 , the dashed line is expected, 
and is a better match to the data. 

It is interesting to check the result as a func- 
tion of position in the detector because (a) some 
reconstruction algorithms are less certain for ver- 
tices close to the PMT wall, and (b) possible 
backgrounds, such as neutrons from nearby cos- 
mic rays [ [jlp^ , or entering events that evade 
the outer detector veto, would accumulate near 
the fiducial boundary. Figure |^a shows R versus 
distance from the PMT wall, where the fiducial 
volume is found at 2 meters. Figure §b shows the 

6 Caveat: the 1MB collaboration used the opposite defini- 
tion, so down-going neutrinos are near cos = — 1 in their 
publications. 
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Figure 1. The zenith angle dependence of R 
for sub-GeV and multi-GeV atmospheric neutrino 
samples from Super-Kamiokande. The dashed 



line shows the expected shape for i/„ 



oscil- 



lation with sin 2 29 = 1 and Am 2 = .005 eV 2 /c 2 



zenith angle dependence of R, dividing the data 
into two approximately equal fiducial volumes: an 
outer volume between 2-5 m from the PMT wall 
and an inner volume greater than 5 m from the 
wall. There is no variation of the result due to 
the fiducial boundary evident in either figure. 

The double-ratio of (jij &)data to (fjt/e)MC is 
useful to illustrate the effect, but it does not indi- 
cate whether /i or e rates (or both) are affected. 
Furthermore, R is not so practical for statistical 
tests. Figure H shows the /j, and e rates separately 
for sub-GeV and multi-GeV (FC+PC) compared 
to Monte Carlo prediction. The solid bands are 
the absolute prediction, where the height of the 
band is equal to the Monte Carlo statistical un- 
certainty. Not shown is the ±20% normaliza- 
tion uncertainty, which is highly correlated bin- 
to-bin, between fj, and e, and between sub-GeV 
and multi-GeV. Even accounting for this uncer- 
tainty, it is apparent that the anomalous R is 
dominated by a deficit of //-like events coming 
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Figure 2. (A) R for the sub-GeV samples 
as a function of distance from the PMT wall 
(Dwall)- (B) R versus zenith angle for two con- 
centric fiducial volumes: 5 > Dwall > 2 meters 
(outer) and Dwall > 5 meters (inner). 



from below (cos 6 < 0). 

The significance of this result can be easily 
evaluated by calculating the up-down asymme- 
try A = (N down - N up )/(N down + N up ) where 
up and down are defined by cosG < —.2 and 
> .2 respectively (Tab.^). Besides other interest- 
ing possibilities |20f| , the distribution of A is nicely 
described by a gaussian variance. For multi-GeV 
events, N up and Ndown should be nearly symmet- 
rical, whereas A for /x-like data (FC+PC) differs 
from expectation by greater than 5a. 

Table 4 



Up-down asymmetry for multi-GeV data. 





S ample 


1 v up 


*V down 


A 


M 


-like data 


102 


195 


0.313±0.055 


e- 


like data 


76 


90 


0.084±0.077 


/' 


-like MC 


1669 


1707 


0.013±0.017 


e- 


like MC 


596 


589 


-0.006±0.029 



The dashed line in Fig]3j represents an oscil- 
lation hypothesis-] of (sin 2 29 = l,Am 2 = .005) 
for disappearance. The overall normalization 
is adjusted upward (thus the e-like rate increases 

7 The values (1, .005) represent a test point; the exact best 
fit location can change when the technique or data sample 
changes since the minimum is fairly flat. 
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Figure 3. The number of /^-like and e-like events 
as a function of zenith angle. The solid his- 
tograms are the Monte Carlo expectation with no 
neutrino oscillation; the thickness represents the 
statistical uncertainty in the Monte Carlo sam- 
ple. The dashed line shows the expected shape 
for — > v T oscillation with sin 2 2f5 = 1 and 
Am 2 = .005eV 2 /c 2 . 



to better match the data, even when v e mixing 
is not considered) while (2) is used to calculate 
the probability of disappearance and reweight 
the Monte Carlo. The v travel distance L is cal- 
culated as a function of energy and flavor based 
on a production height model[^3|. The oscillation 
hypothesis provides a reasonable fit to the data, 
certainly better than the null hypothesis. 

The exact details of fitting the data to esti- 
mate possible mixing parameters are still being 
evaluated. Using a method similar to that used 
by Kamiokande fl5|l , x 2 terms are formed between 
the data and Monte Carlo prediction [modified by 
P U u'(sm 2 29, Am 2 )], binned in zenith angle, en- 
ergy, and flavor (values of R are not used directly 
in the fit). The normalization, N^/Ne ratio and 
systematic terms are allowed to adjust and con- 
tribute to the x 2 - For the Super-Kamiokande 
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— 90% C.L. 
99% C.L. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

sin 2 29 

Figure 4. Confidence intervals for sin 2 29 and 
Am 2 based on a x 2 fit to Super-Kamiokande at- 
mospheric neutrino data binned by lepton identi- 
fication, lepton momentum and cos0. The solid 
line is 90% CL, the dashed line is 99% CL. 



data, the minimum \ 2 is found to be rather likely, 
~30% depending on the details. From the fit, a 
confidence interval is drawn based on Xmin + 4.6 
(90% CL) such as shown in Fig. § 

The preferred interval for sin 2 29 is found near 1 
because the upward zenith bin with (L) ~ 10, 000 
km, presumably has averaged over several oscil- 
lation lengths and P vv < is one half. In a scenario 
where more than two v flavors are mixing, the av- 
erage value can be less than one half. The Am 2 
range is determined by the shape of the zenith 
angle, also considering dependence on E v . This 
90% confidence interval from Super-Kamiokande 
prefers a lower range in Am 2 than that found by 
the Kamiokande collaboration[[l5|, which had a 
minimum Am 2 of .005. 

6. Upward- Going Muons 

The above discussion covered neutrino interac- 
tions in the fiducial volume of the detector. The 
other class of atmospheric neutrino event studied 
is that of interactions in the rock around the 



detector, where the final state muon enters the 
sensitive region of the detector. To separate these 
from ordinary cosmic ray muons, the muon must 
be upward-going, or come from the direction of 
a known thick overburden. The parent neutrino 
energy is 10-1000 GeV, significantly higher than 
for contained events. 

There are several current measurements of the 
total rate of upward-going muons. In addition 
to the water Cherenkov detectors described else- 
where in this paper, MACRO and Baksan are 
large area scintillator detectors that distinguish 
upward-going muons by time-of-flight; T. Mon- 
taruli (MACRO) and S. Mikheyev (Baksan) have 
presented updated results at this conference. The 
measured and expected event rates are compared 
in Table ||; the experiments are listed in order 
of increasing absorber depth for directly vertical 
muons. The uncertainty in the measured num- 
ber of events is statistical only; some experiments 
have estimated that the uncertainty due to exper- 
imental systematics could be as large as 8%. The 
uncertainty in the expected number of events is a 
common 15-20%, due mostly to the uncertainty 
in the absolute flux. Considering this, in no case 
is a significant deficit of muons measured, but the 
measurements are generally low compared to ex- 
pectation. 

Table 5 

Summary of through-going upward-// totals. 

Number of Events 



Experiment Measured Expected 



MACRO ]24| 


350 


± 


19 


472 


Baksan |25|] 


558 


± 


24 


557 


Kamiokande p6| 


373 


± 


19 


414 


IMB|27| 


539 


± 


23 


550 or 625 


Super-Kpf 


410 


± 


20 


445 



There are two other approaches that probe neu- 
trino oscillations |2j| using upward-going muons. 
1MB measured the ratio of stopping upward 
muons {{E v ) ~ 10 GeV) to through-going upward 
muons ((E v ) ~ 100 GeV) to be 0.16 ± .02. Al- 
though much of the flux uncertainty cancels out, 
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Upward Through Going Muon Flux 



the usual analytic integration using deep inelastic 
scattering and parton distribution functions must 
be handled with care l^JjOf , especially at low 
energy After these considerations, the pre- 
dicted rate was .14 or .18 depending on the flux 
model|27|]. Based on the agreement of data with 
prediction, a small excluded region in sin 2 29, 
Am 2 was drawn around Am 2 = 10~ 3 eV 2 /c 2 , 
where the strongest deviation would have been 
found; this happens to be in conflict with the 
region favored by the Super-K contained vertex 
data. Super-Kamiokande will also measure the 
stopping ratio, with the advantage of a very thick 
detector that stops a large number of upward- 
going muons; the analysis is in progress. 

The second approach uses the shape of the 
zenith distribution, which may be distorted as 
the baseline varies from 500 km at the horizon to 
12,000 km at the nadir. Because the energy spec- 
trum of the parent neutrinos is broad, ~ 10 — 1000 
GeV, the change in shape is gradual, with some 
steepening at the horizon as the probability de- 
creases for high energy to oscillate. Figure || 
shows the Super-K measurement of the flux ver- 
sus zenith angle compared to expectation; recall 
that the normalization of the prediction is uncer- 
tain to ±20%. When the normalization of the 
Monte Carlo is decreased by a factor of a = 0.83 
to match the data, the \ 2 is 12.7; alternatively, 
when the normalization is increased by a factor 
of a = 1.2, and disappearance oscillations are 
applied with sin 2 28 = 1, Am 2 = .005, the fit is 
somewhat better, with a x 2 of 8.3. 

The other experiments listed in Table |B| can 
make this comparison as well. The scintillator de- 
tectors MACRO and Baksan unfortunately have 
reduced acceptance at the horizontal, so those 
bins require significant geometric correction. Cu- 
riously, even the well-measured upward bins of 
those two experiments are not smooth and suffer 
from poor agreement in shape with: (a) no os- 
cillation, (b) any 2-flavor oscillation parameters^, 
and (c) each other. Both experimental groups 
have made extensive checks for a systematic er- 



8 A recent preprint considers that matter oscillation with 
a sterile neutrino crossing the Earth's core may modu- 
late the prediction with features similar to the MACRO 
data||]. 
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Figure 5. The through-going upward-going fx 
flux as a function of zenith angle as measured 
by Super-Kamiokande. The data points are com- 
pared against expectation (solid histogram) and 
the same expectation with normalization a = 1.2 
and ^-mixing sin 2 28 — 1, Am 2 = .005. 



ror, but have found no cause |25|j32(| . 

7. Conclusions 

Currently, the evidence for an anomalous ra- 
tio of atmospheric neutrino flavors is inconsistent 
across experiments. Of course, prior results re- 
main intact, and either support or disagree with 
the anomaly. But the latest results support that a 
significant discrepancy exists between experiment 
and prediction. The measurement by Soudan-2 
shows that the anomaly is not specific to water 
Cherenkov detectors. A zenith angle measure- 
ment of the Soudan-2 data could be very inter- 
esting. Significant new information is taken from 
the high statistics Super-Kamiokande data: the 
anomaly is strongly confirmed in R and the zenith 
angle dependence of R. The shape of the zenith 
angle dependence is very suggestive of neutrino 
oscillations. The multi-GeV /i-like rate as a func- 
tion of zenith angle indicates that disappear- 
ance is favored over — > v e oscillation. Even 
though 1.1 years of Super-Kamiokande running 
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has surpassed the prior world statistics, more live- 
time will be welcome and allow finer subdivision 
of the data for cross-checks and estimation of pos- 
sible mixing parameters. 

Still, to claim the atmospheric v anomaly is 
caused by neutrino mixing requires confirma- 
tion with as many other signatures as possi- 
ble. Upward-going muons currently show a mix 
of results, none obviously refuting or confirm- 
ing. The upward /i flux measurement of Super- 
Kamiokande is consistent with the estimated pa- 
rameters from the contained vertex sample, but 
the result is not conclusive. Further analysis of 
stopping muons and the angular distribution is 
awaited. 

In the end, the best confirmation should come 
from long-baseline neutrino beams from accel- 
erators. However, according to the preliminary 
Super-K results, the experiments under construc- 
tion may have to address lower values of Am 2 
than when they were conceived. 

8. Acknowledgements 

The author thanks the conference organizers 
and the participants listed in this paper. He 
is especially appreciative of his collaborators on 
Super-Kamiokande for preparing the latest re- 
sults, which are presented here on their behalf. 
He gratefully acknowledges financial support by 
the U.S. Department of Energy. 

REFERENCES 

1. T.K. Gaisser et ai., Phys. Rev. D54 (1996) 
5578. 

2. M. Honda et ai., Phys. Lett. B248 (1990) 193; 
M. Honda et ai, ibid.T)52 (1995) 4985. 

3. G. Barr et ai., Phys. Rev. D39(1989) 3532; 
V. Agrawal et ai., ibid.D53(1996) 1314; 
T.K. Gaisser and T. Stanev, Proc. 24th Int. 
Cosmic Ray Conf.(Rome) Vol.1 (1995) 694. 

4. M. Circella et ai., Proc. 25th Int. Cosmic Ray 
Conf.(Durban) Vol.7 (1997) 117. 

5. G. Barbiellini et ai., Proc. 25th Int. Cosmic 
Ray Conf.(Durban) Vol.6 (1997) 317. 

6. G. Tarle et ai., Proc. 25th Int. Cosmic Ray 
Conf.(Durban) Vol.6 (1997) 321. 



7. T. Stanev, These proceedings. 

8. Y. Fukuda et ai., sub mitted to Phys. L ett., 
ICRR-Report-41 1-98-7, |hep-ex/9803006| . 

9. D. Casper et ai., Phys. Rev. Lett. 66 (1991) 
2561; R. Becker-Szendy, et ai., Phys. Rev. 
D46 (1992) 3720. 

10. K.S. Hirata et ai., Phys. Lett. B205 (1988) 
416; K.S. Hirata et ai., ibid.B280 (1992) 146. 

11. W.W.M. Allison et ai., Phys. Lett. B391 
(1997) 491 

12. T. Kafka, These proceedings. 

13. K. Daum et ai., Z. Phys. C66 (1995) 417. 

14. M. Aglietta et ai., Europhys. Lett. 8 (1989) 
611. 

15. Y. Fukuda et ai., Phys. Lett. B335 (1994) 
237. 

16. R. Clark et ai., Phys. Rev. Lett. 79 (1997) 
345. 

17. J. Engel et ai., Phys. Rev. D48 (1993) 3048. 

18. K. Inoue, These proceedings. 

19. S. Kasuga et ai., Phys. Lett. B374 (1996) 238 

20. J. Flanagan et ai., UH-511-880-97, |he£- 



24. 
25. 
26. 



27. 

28. 

29. 

30. 
31. 

32. 
33. 



ph/9709438 



21. O.G. Ryaznskaya, JETP Lett. 61 (1995) 237. 

22. Y. Fukuda et ai., Phys. Lett. B388 (1996) 
397. 

23. T.K. Gaiser and T. Stanev, BRI-97-28, |astro 



ph/9708146 



T. Montaruli, These proceedings. 
S. Mikheyev, These proceedings. 
Y. Totsuka, 28th Int. Symposium on Lepton 
Photon Interactions, Hamburg, 1997; Kamio- 
kande collaboration, To be published. 
R. Becker-Szendy et ai., Nucl. Phys. B38 
(Proc. Supp.) (1995) 331. 
Super-Kamiokande collaboration, To be pub- 
lished. 

P. Lipari and M . Lusignoli, ROME-1190-97, 
hep-ph/9712278| . 

W. Frati et ai., Phys. Rev. D48 (1993) 1140. 
P. Lipari et ai., Phys. Rev. Lett. 74 (1995) 
4384. 

MACRO collaboration, In preparation. 



QY. Liu and AYu. Smirnov, IC-97-211, |hcp- 



ph/9712493 



